A=[1 0 1; 0 1 2; 1 1 1];
x=[1;4;7];
inv(A)*x
